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We investigate the link between information and thermodynamics embodied by Landauer’s principle in the 
open dynamics of a multipartite quantum system. Such irreversible dynamics is described in terms of a col- 
lisional model with a finite temperature reservoir. We demonstrate that Landauer’s principle holds, for such a 
configuration, in a form that involves the fiow of heat dissipated into the environment and the rate of change 
of the entropy of the system. Quite remarkably, such a principle for heat and entropy power can be explicitly 
linked to the rate of creation of correlations among the elements of the multipartite system and, in turn, the 
non-Markovian nature of their reduced evolution. Such features are illustrated in two exemplary cases. 


Logical irreversibility and heat dissipation are linked 
through the relation provided in 1961 by Landauer HI: the 
erasure of information on the state of a system is concomi¬ 
tant with the dissipation of heat into the surroundings. In turn, 
such heat is lower-bounded by a change in the information- 
theoretic entropy of the system. Landauer’s principle, which 
has been recently tested m , was the building block of remark¬ 
able advances in our understanding of thermodynamics and its 
links with logical irreversibility and information theory O. 

Recent progress in the quantum approach to non¬ 
equilibrium statistical mechanics Elia has made the tracking 
of quantities such as heat, work, and entropy possible in an ex¬ 
perimentally viable way. In turn this has enabled the test-bed 
demonstration of the link between information and energy in 
quantum systems subjected to elementary quantum computa¬ 
tion protocols m. The relation between information and the 
thermodynamic costs of quantum operations has been exten¬ 
sively addressed in the recent past. Noticeable examples in¬ 
clude Refs. Q, (cf. in for a recent overview). Techniques of 
quantum statistical mechanics have been used to prove that a 
finite-size environment can provide tighter bounds to the heat 
generated in an erasure process Enoi. The validity of Lan¬ 
dauer’s principle has been verified for a quantum harmonic 
oscillator strongly coupled to a bath of bosonic modes ifTTI . 
However, such conceptual and experimental progress has not 
yet resulted in a satisfactory microscopic quantum framework 
able to account for the emergence of Landauer’s principle. 

In this Letter we provide a derivation of a Landauer-like 
principle that addresses the heat and entropy fluxes implied 
in an erasure process described by a collision-based picture 
of open-system dynamics (121 . Such a microscopic formu¬ 
lation has been used to examine the process of thermaliza- 
tion of a quantum system in contact with a non-zero tem¬ 
perature bath and to investigate the link between with non- 
Markovianity and quantum correlations iHMa. Here, we 
bring together these two perspectives: We analyse the process 
of information-to-energy conversion in collisional models that 
make use of general, non-restrictive assumptions and effective 
mechanisms to describe the heat dissipation into an environ¬ 
ment. Our study thus aims at bridging the gap highlighted 
above and contributes to the ongoing research for a satisfac¬ 


tory foundation of Landauer’s principle. The advantages of 
our approach are manyfold. First, modelling the irreversible 
dynamics of a system by a collisional model where, at each 
step, the environment is reset in a thermal state enables the 
formulation of Landauer’s principle in terms of (heat and en¬ 
tropy) powers. Second, we show that the amount of heat dissi¬ 
pated by a multipartite quantum system is lower-bounded by 
the correlations established among the elements of the sys¬ 
tem itself. Finally, the flexibility of the picture adopted here 
allows us to unveil interesting links between the violation of 
Landauer’s principle and the emergence of memory-bearing 
effects leading to non-Markovian dynamics. 

Erasure through thermalization- We focus on the thermaliza- 
tion process of a system S in contact with an environment ^ 
endowed with a large number of degrees of freedom, which 
can thus be treated as a bath ca. In our model R consists of a 
collection of N identical non interacting elements hence 
after dubbed subenvironments, each of which is assumed to be 
in a thermal state. The bath is therefore in the product thermal 
state r|= fl* with r|* = the 

free Hamiltonian of the n^^ subenvironment, and P the inverse 
temperature. We call Hy the free Hamiltonian of the system 
and, with no loss of generality, we assume for any 

n. The system interacts with the environment via a sequence 
of pairwise collisions with individual subenvironments. The 
assumption of a big reservoir implies that the system never 
interacts twice with the same subenvironment, therefore at 
each collision the state of the subenvironment is refreshed and 
the system always interacts with subenvironments in a ther¬ 
mal state. Each collision is described by a unitary operator 
U = where g is a coupling constant, x is the collision 

time and y=Y.k^k^^k is the (dimensionless) S-%^n interac¬ 
tion Hamiltonian [S^ (Rk) are Hermitian operators acting on 
S (^n)]- The information initially encoded in p is gradually 
diluted into After the (^-1-1)* collisions, we find the states 
of the system and r|„+i of the environment 

Pn+i ~ Tr^[0^ (8) = d^[p^], (1) 

fln+i ~ Tr^ [U = A^[t|^], (2) 

where d> (A„) is a completely positive trace-preserving 


2 


(CPTP) map on S {%) (HI. Unlike map does depend 
on n through p„. A formal definition of mean heat flux and 
work is obtained considering infinitesimal changes of internal 
energy as dTr \juH] = Tr [judH] +tr[5^^], where ju is the den¬ 
sity matrix of a system with Hamiltonian H (Tfill . The first 
term is identified with the average work done on/by the sys¬ 
tem; the second, which disappears for a unitary evolution, is 
a heat fiux. Such arguments have been more recently reprised 
in Refs. (13. In this sense, based on Eqs. 0 and the 
corresponding energy variation of the system AE^+i and heat 
exchange with reservoir AQ^+i, are given by 


A£„+i =Tr [H5(4>-I)[p„]], AQn+i =Tr [H 3 ^(A„-I)[ii]] . 

We next assume (R)^)ri^=Tr^[Ry^r|;^]=0 (this is possible by 
moving into an interaction representation with respect to a 
suitably redefined local Hamiltonian of S and without affect¬ 
ing our results). Assuming that each collision lasts a short 
time, we can expand (<I> — I) in series up to second order 
in y = gx to find Ap„+i = (<!>- I)p;^ = ^pn with ^ = 
Tr^[y(p„ (g)r|)y - {y^, (g)r|}/2]. This gives us 

AE„+, = gh^'£i^kRMSASj-H^kSj,Hs})p„, (3) 

kj 

AQn+l = g^^'^Y.(^kSj)p„{RkH^Rj-H^kRj,HR})r^. (4) 
kj 

In the limit x 0 and 1, the quantity t = nx becomes 
a continuous variable, ^ p(r) and (p„+i — p„)/x ^ p, so 
that we get the master equation (ME) p = ^[p(t)]. Likewise, 
Eqs. 0 and 0 become 

E = yY.{RkRjU{SkHsSj-\{SkSj,Hs})p , ( 5 ) 

kJ 

Q = y'£{SkSj)p{M^Rj-HRkRj,HR})r^, (6) 

kJ 

where we have defined the rate y = g^x, achieved by taking 
both g ^ oo and x ^ 0 so that g^x is a constant. By doing so, 
we would guarantee the system-environment effects to per¬ 
sist (differently from what is obtained by taking x ^ 0). As 
discussed in Refs. (ISllMl, with the assumption on {Rk)y] in¬ 
voked above, such a continuous-time limit is fully consistent 
and legitimate. By construction, p(r) evolves according to a 
CPTP dynamical map (121 [ID. We now focus on energy- 
conserving S-^ interactions and assume [U,{Hs-\-H^)]=0. 
Hence, Q = —E and the stationary state of the system is the 
Gibbs state at the same initial temperature of the bath, i.e. 
p^^=e~^^s In this case, the stationary state is 
such that E = Q = 0, which is not verified by non-energy- 
conserving models. Is such cases, the steady state corresponds 
to a constant fiux of heat into/out of the system. 

The relative entropy ^(plp^^) between the state at time t 
and the stationary one obeys the relation 5(p|p^^) = —5(p) + 
(m. Due to the non-increasing nature of the relative en¬ 
tropy under CPT maps (20l , we obtain 

pe(0>5(p), (7) 


where S{p) = —S{p). This embodies an open-system formu¬ 
lation of Landauer’s principle for the heat and entropy fluxes. 
We can thus lower-bound the flow of heat dissipated into the 
environment with the rate of change of a key information the¬ 
oretical quantity. Remarkably, our formulation enables the as¬ 
sessment of the nature of the open-system dynamics of 5. The 
power of our formalism is thus twofold: on one hand, it allows 
for a time-resolved analysis of the erasure-by-thermalization 
process. On the other hand, it elucidates the role of correla¬ 
tions in the efficiency of information-erasure processes. This 
is illustrated by two examples of irreversible collisional dy¬ 
namics. We first analyse a cascaded system where the suben¬ 
vironments interact in sequence with the elements of a bipar¬ 
tite system, chosen here for easiness of presentation. In this 
scenario, the environment is “recycled” to erase the informa¬ 
tion in the two subsystems before being damped. We show 
that the efficiency of information erasure is limited by the 
build-up of intersystem correlations. In the second example, 
which we dub indirect, we consider a system interacting with 
an ancilla whose state is then erased by the environment. In 
this case, we show the appearance of local violation of Lan¬ 
dauer’s Principle due to the onset of memory effects. 
Information erasure in cascaded systems.- Eor now, we con¬ 
sider the system S as bipartite and consisting of two identical 
subsystems 5 a, 5 - Each subenvironment collides with 5 a 
first and Sb later. Each collision is modelled as a unitary pro¬ 
cess with associated evolution operator Ux=e~^^^^^, (X=A^B) 
and generator Vx='Lk Sxk^Rk (here Sxk acts on subsystem X). 
A “step” consists of an Sa-%i collision followed by an 5s- 
one. Hence, the joint 5-^ state at step ^+1 is given by 
On^l=UBUApn^ulul (pn is the joint 5 a- 5s state). Eqs. ^ 
and ^ thus still hold with U = UbUa- 

Through a procedure analogous to the one outlined earlier, 
it can be shown that in the continuous-time limit the bipartite 
system 5 is governed by the ME da 

P= X 2A[p]-f ^2 )a5[p], (8) 

X=A,S 


where superoperators Lx and Dab are defined by 

^x[p] (9) 

kj 

'^AB[p]=y'^{RjRk)^ [‘^Ay^P^^^Sj] [‘^Sj, P^^A)^] -(10) 

kj 


Subsystems 5a and Sb jointly undergo a Markovian dynamics. 
This is true also for the dynamics of 5a, given that it collides 
with “fresh” subenvironments every time ED. However, in 
general, 5s is subjected to non-Markovian evolution. 

As before, given that no work is performed on the system, 
we identify the variation rate of the internal energy of the sys¬ 
tem E with AE^+i = Tr [^j(<I>—I) [p„]]. Using Eq. (15) and 
following an approach similar to the one presented above, it 
is possible to write explicitly the expression for the heat trans¬ 
ferred at each step of the collision-based process, as shown 
in the Appendix. By focusing again on energy-conserving 
system-environment interactions, the identity Q = —E holds 
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FIG. 1. (Color online) (a) Sketch of the recycled environment setup: after colliding with Sa^ subenvironment (prepared in is used to 
erase the state of Sb- By tracing out the subenvironments, a Markovian ME is achieved for 5. (b) Heat flux and entropy for a system composed 
of a single qubit for ^ = 0.9 and ^ = —0.9 [inset], (c) Heat and entropy flux for the cascade model. We have set y/co= 1 with S initially 
prepared in |tt) at ^ = 0.9 and ^ = —0.9 [inset], (d) Mutual information /(A : B) between the subsystems in the same case of panels (b) and 
(c) for ^=0.9 (solid brown line) and —0.9 (dashed green line). Heat fluxes are expressed in unit of coy. 


and the system thermalizes at the same temperature of the en¬ 
vironment. Accordingly, in light of the absence of any direct 
SaSb coupling, the system asymptotically reaches the sta¬ 
tionary state (8) p^^ = / (ZaZb) (here 

Zx = ]), i.e. the a factorised state of locally thermal 

states at the same temperature. Indeed, despite the Sb evolu¬ 
tion depends on Sa (thus determining its non-Markovian char¬ 
acter), Sa thermalizes with ^ in a fully Markovian fashion. 
Thereby, asymptotically, the Sa-%i collisions do not change 
the state. Thus, in the same limit, Sb will also collide with 
subenvironments that are still in state r|, forcing it to thermal- 
ize at the temperature of the environment. We remark that the 
bound stated in Eq. (|7]) clearly holds in the present situation, 
provided that Q and ^(p) refer to system 5. 

In order to show that recycling the environment affects 
the erasure efficiency, we use the definition of quantum mu¬ 
tual information I{Sa '• Sb) = S{Psa)~^^{PSb) “‘^(p) con¬ 
ditional quantum entropy =‘5(pM(PAt) C21, to get 

-5(p) = -25(p5^) + and, in turn 

pe > 25(p5j +5^,15, +/(5 a:5b). (11) 

This shows that the erasure efficiency in cascaded systems de¬ 
pends on the rate at which correlations between Sa and Sb are 
established. Indeed, all the quantities in Eq. ( pT] ) are null only 
for product states. Incidentally, the rate 25(pj^) is that of two 
identical state erased by independent reservoirs. 

Our arguments are straightforwardly extended to the multi¬ 
partite case. Using the correlation information St — 

Si...N 1^ . which extends mutual information to an A—partite 
system, Eq. ( pTj ) is extended to the erasure of N copies as 

Pe>;V5(pi) + (lV-l)5l|2,.,;V-i:L2‘5^|l..+^V. (12) 

We now illustrate the results described until now with a spe¬ 
cific example. We consider an all-qubit configuration where 
Hs=M^)Lx=a,b^Xz and with {da} (/ = 

x,y,z) the Pauli operators of qubit ^ = Each subenvi¬ 

ronment is prepared in the thermal state r|=[( 1 — ^)/2] |t) (t I + 
[(l+^)/2]||)(||, where are the eigenstates of 

with ^=tanh(Pco/2). The Sx-%n interaction Hamiltonian is of 
the isotropic XX form, i.e., Vx='Li=x,y^xi ^ In the Ap¬ 
pendix we present the form taken by Eq. ( p~5] ) in this case and 


show that the evolution of Sa is Markovian and not causally 
related to Sb- Instead, the Sb dynamics depend crucially on 
the state of Sa, thus being steered by previous subsystem- 
environment collisions. The behaviour of P2 is compared to 
5(p) in Eig. for two different temperatures and each sub¬ 
system initially prepared in |t). To pinpoint the role played 
by the SaSb correlations in setting the non-monotonic tem¬ 
poral trend of the heat fiow [cf. Eig.[^(c)], we have studied 
I {A : B). Eig.[^(d) shows that the establishment of such cor¬ 
relations indeed acts as a precursor of the non-monotonicity 
of the heat exchanged with the environment. In the Appendix, 
we support these conclusions with the study of a multipartite 
system. 

Indirect Erasure- We now analyse the dynamics of a bipartite 
system 5, in which the evolution of the two subsystems and 
their coupling with the subenvironments are treated on dif¬ 
ferent footings: we include a mutual interaction between Sa 
and Sb while we assume that only Sb interacts with This 
amounts to assume that the information contained in Sa is first 
transferred to Sb and then damped into a clear example of 
indirect information erasure. In practice we assume that be¬ 
tween the Sb-%ji and SB-%^n+\ collisions, the two subsystems 
Sa and55 interact. Eqs. 0 and 0 are still valid provided that 
tj UbUab, where Ub (describing an Sb-^ collision) is the 
same as in the last Section, while Uab = regulates the 

joint evolution of the SaSb system. To go to continuous-time, 
we expand each unitary operator to lowest-order in x. Unlike 
the previous case, we should now take Ub — I—, 
as before, and Uab — I—i'iHab- We thus find 

P= — —(13) 

kj 

This is a Lindblad ME describing Markovian dynamics of 5. 
The dynamics of Sa, however, is in general non-Markovian. 

We now illustrate the indirect erasure model with a specific 
example where we set no constraint on the dimension of the 
Hilbert space of Sb- In the Appendix, we derive the form of 
Eq. ( p^ when Sb is a quantum harmonic oscillator resonantly 
coupled to a two-level subsystem Sa via the interaction Hamil¬ 
tonian HAB=(^{q^^P^+^Az/'^)-^J{4^Ax+P^Ay), with q and p 
the position and momentum quadrature operators of oscillator 
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FIG. 2. (Color online) (a) Sketch of the indirect erasure protocol: the bipartite system S consists of subsystems Sa and Sb- After Sb interacts 
with the subenvironment (prepared in rj*), it collides with Sa and is then directed to element Sb thus bridges the erasure process 

undergone by Sa^ which experiences a non-Markovian evolution, (b)-(e): Total heat flow pg against the entropy change rate S for a qubit-qubit 
(b,c) and qubit-harmonic oscillator (d,e) configurations. We have set //co = 0.1 and 7^/00 = 0.01, taking Sb at the same temperature of the 
environment, P=10 (b,d) and P=0.5 (c,e). A study against the value of J and jg is reported in the Appendix. 


Sb, respectively. An analogous model holds for the Sb — 
interaction, each subenvironment being a qubit. The bosonic 
nature of Sb allows us to explore the effects that multiple ex¬ 
citations in Sb have on both the heat flow and rate of entropy 
change. Indeed, as illustrated in Figs. (d) and (e), for an 
initially thermal harmonic oscillator at the same temperature 
of the bath, the differences with a flnite-dimensional subsys¬ 
tem Sb are striking, and dependent on p. At low temperature, 
due to the excitation-conserving nature of the Hamiltonian, 
the dynamics of S is in fact identical to that of a two-qubit 
system, as only one excitation at most can be exchanged with 
This is not the case when the temperature is raised, as ex¬ 
citations can be pumped into the harmonic oscillator at each 
collision [cf. Fig. [^(e)]. In the Appendix we study the be¬ 
havior of the key quantities in our analysis against J and y^. 
Interesting considerations on the nature of the dynamics of 
Sa are in order. Indeed, by tracing out the degrees of free¬ 
dom of Sb, we obtain a non-Markovian evolution due to the 
flowback of information into Sa- Correspondingly, Eq. 0 is 
violated very quickly, as shown in Fig. which again high¬ 
lights the effects that the establishment of intra-system cor¬ 
relations have on the efficiency of the information-to-energy 
conversion processes. Similar conclusions hold when relaxing 
the key assumption of non-interacting subenvironments. The 
break-down of such condition would entail the establishment 
of memory effects within the environment and the occurrence 
of non-Markovianity in the SaSb dynamics. The violation of 


(a) (b) 



FIG. 3. (Color online) We plot P2 a( 0 and Sa for the qubit-oscillator 
configuration in the indirect-erasure case. In panel (a) [(b)] we have 
used P = 10 [P = 0.5] and the parameters stated in Fig.|^(d) and (e). 


our formulation of Landauer’s principle could thus be used to 
infer the nature of the bath itself. 

Conclusions. We have investigated the information-to-energy 
conversion processes in terms of collision-based models to de¬ 
scribe the system-environment interaction, with the aim to go 
towards the a microscopic formulation of Landauer’s princi¬ 
ple. In our approach we focused on the rate of entropy change 
and dissipated-heat flux involved in the overall erasure pro¬ 
cess. By doing so, we have been able to link the efficiency of 
the erasure process to the intra-system correlations arising in 
the open dynamics of a multipartite system. 
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APPENDIX 

Heat flux in the qubit ‘recycling environment’ model 

As outlined in the main text, no work is performed on the 
system, so we identify the variation rate of the S internal en¬ 
ergy E with 


Tr 

A£'„+i = lim — 

n^oo 

T^O 

Remembering the master equation 
P= 51 -^A-[P]+ ®Ab[P], 

X=A,B 

Ax[p]=yY.{RkRj)T\{SxjPSxk-^{SxkSxj,p}), 

kJ 

‘DAB[p]=y^{RjRk)^ [5a<:P,>5b;] +{RkRj)^ [^By'jP^A*:] -(15) 
kj 


[HA<&-I)[P„]] 


(14) 
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this reads 


E = Tr 

i T 2:x[p]+®ab[p] 

) E H,, 


_ \X=A,B ^ 

/X=A,B 


with the free Hamiltonian of subsystem X=A,Follow¬ 
ing the recipe presented in the main text, it is possible to write 
explicitly the expression for the heat transferred at each step 
of the collision-based process, this is given by the quantity 
A2„=Tr[^^(r|„ — p)]. When taking the continuous limit of 
oo and x ^ 0, we find that the rate of change of the heat 
exchanged with the environment reads 


E = y 


Y,^kj{{SAkSAj)p + {SBkSBj)p) ^'^^^{?>kj{§AkSBj)p) 
kj 


with akj={RkHRRj-^{RkRj,HR})r^ and ^kj={Rj[Hii,Rk\)r\- 


Example I: Master equation for the qubit ‘recycling 
environment’ model 


l3Q{t)/N §N{t)/N 



FIG. 4. (Color online) Recycling environment configuration for dif¬ 
ferent numbers of subsystems N. We have used y/co = 1 with ^ = 0.9. 
All the subsystems Sn are prepared in |t). All the quantities, i.e. heat 
flux pg, entropy variation 5, mutual information and radiation rate 
(/+/“), have been rescaled with N for better visibility 


As described in the main text each collision is generated 
by a Hamiltonian of the form Hx = 'Lk^f^^k- Here, = 
{af and R = We assume the sub environ¬ 

ment ^ initially in a thermal state at temperature p, defin- 

_P 

ing ^=tanh(p0)/2), this means 

Tv[e~^2^z ] 

Introducing the Pauli ladder operators a^=(d^±/d^)/2 for 
Q=SA,B^%.n and after straightforward calculation. The master 
equation in Eq. (14) of the main Letter takes the form 



p' =L^[p*]+L^[p‘] + ‘D^^[p‘] 

L^[p‘] =r+L[a-J(pO + r-L[a+J(p' 

‘D^^[p‘] =r 

r 



A 

- 


"a) 



- 




with p^ the density matrix of the overall system, y = g^x, 
= 2y(l ±^) and L[o\{p)=6pd^ — ^{66^,p} for a generic 
operator d. Needless to say, the reduced dynamics of Sa is 
not causally linked to that of B due to the temporal order¬ 
ing of the interactions between the sub-systems and a given 
sub-environment. The opposite is, quite obviously, not true. 
Explicitly 

Pa = r+L[a-J p^^+r-L[a+J p^^, 

Pa (17) 

where the expectation value in the second line of Eq. ([Tt]. is 
evaluated over the reduced state of sub-system 5 a- The ex¬ 
pression for the heat flow from the overall system 5 can con¬ 
veniently be given in terms of the elements = (/,7|p^|k,/) 


(/, 7 , k, I =t, D of the two-particle density matrix. Explicitly 


Q{t) = 4y(0 1 + pti + pti) + (p- 


ft ^ 





(18) 


with lit), II)} the eigenstates of d^. Eq. ( p^ shows the direct 
dependence of the heat fiow on the coherences of the 5 state. 
Similarly, 


QA{t)—2y(0 (l-l-^)(pjj+pj|)-(l-^)(p|j+p||) 
Qb{ t) =2yo) (1 +^) (p||+p|j) - (1 -^) (p||+p||)+ 

2 ^(pIt + PtI) 


(19) 


Extending to a multipartite system 5 composed by N qubits, 
the master equation becomes 


p' = fx'^[p'] + f ®'^'[py (20) 

i j>k 

In Eig.Q we report the heat flux and entropy variation be¬ 
haviours for recycling environment in cases up to 5 qubits. 
Eor = 1 this is the heat fiow for a single qubit coupled to 
the environment with no recycling. Looking at the heat fiux 
(Eig|^a)), the initial growth is always greater in magnitude 
and lasts for longer as N increases. This can be explained as 
a cooperative effect between the N qubits which gets stronger 
as N increases. Indeed, the master equation ( [20| can be re¬ 
arranged, writing a non-local coherent part due to a coupling 
mediated by the environment and collective decay channels 
with jump operators = 'Lk^% EH. This effective interac¬ 
tion characterises the initial behaviour of the entropy variation 
and mutual information. 
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Example II: Indirect erasure via a qubit-harmonic oscillator 
system 


We now address the case of indirect erasure achieved 
through the use of a harmonic oscillator bridging the 
interaction of a qubit sub-system Sa and the multipar¬ 
tite environment The interactions between the sys¬ 

tem qubit and the oscillator are described by resonant 
excitation-exchange interactions of the Jaynes-Cummings 
form. Specifically, we call b and P the bosonic annihila¬ 
tion and creation operators of the harmonic oscillator Sb and 
Uab = the time-evolution operator generated by the 

Hamiltonian with 

q={b^b'^)/ and p=i{b^—b) / ^/2 the position and momen¬ 
tum quadrature operators of the oscillator, respectively, and 
b{b^) the corresponding annihilation (creation) operator. This 
accounts for the intra-system collisions. Similarly, the inter¬ 
action between the oscillator and the element of the envi¬ 
ronment is described by the unitary operator Ubr = 
with Hbr given by Hab with the replacements J ^ g and 
^ G^^. By putting in place the procedure for deriving 
a continuous-time master equation, we find that the dynamics 
of system 5 is given by 

( 21 ) 

The heat flow associated with the open-system nature of the 
evolution of 5 can then be expressed as functions of density 
matrix elements ( 7 ,^ and n,m eZ) 

with lit), II)} the eigenstates of G^^ and {|^)} the Fock-state 
basis for Sb- We have 


k=^^ln=0 


( 22 ) 


Heat flux and entropy variations versus coupling constants 



FIG. 5. (Color online) Total heat flow pg (red) and the entropy 
change rate S (blue, dashed), for various values of coupling con¬ 
stant J [in units of co] for the qubit-qubit conflguration. We have 
set y^/co=0.1, taking Sb initially at the same temperature of the en¬ 
vironment, P = 10. 



FIG. 6. (Color online) Total heat flow pg (red) and the entropy 
change rate S (blue, dashed), for various values of rate y [in units 
of co] for a qubit-qubit conflguration. We have set //co=0.5 taking 
Sb initially at the same temperature of the environment, P = 10. 


Here we report the behavior of the heat flux and the entropy 
variations against the coupling rate J between Sa and Sb in the 
indirect erasure protocol, for cusing on the qubit case for defi¬ 
niteness. Fig.|^shows the changes in the temporal behavior of 
both pg and Sab when / = 0.1,0.2,0.4 and 0.8. When investi¬ 
gated within a fixed-size temporal window, the increase of the 
qubit-ancilla coupling induces the appearance of an increas¬ 
ing number of oscillations both in the exchanged heat and the 
entropy changes. This is a clear result of the stronger coherent 
coupling between Sa and Sb, which induces more substantial 
changes in the state of the former. 

On the other hand, an increase in the effective damp¬ 
ing rate y^ resulting from a ramped-up value of the ancilla- 
environment coupling rate g results in a quicker equilibration 
of the state of Sa- the system sees a fasted dissipative dynam¬ 
ics and adjusts to it very quickly. This is shown in Fig.[^ 
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